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We investigate the propagation of light with ultra low group velocity in a Bose-Einstein condensate 
where the phase is not uniform. The light is shown to couple strongly to the phase gradient of the 
condensate. The interaction between the light and the condensate enables us to perform a phase 
imprinting where the phase of the condensate is imprinted on the light. We illustrate the effect by 
showing how one can measure the fluctuating phase in an elongated quasicondensate. 
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Superfiuids such as Bose-Einstein condensates of alkali 
atoms [0 are systems which show a conceptual simplic- 
ity yet provide all the intriguing qualities of an interact- 
ing many-body system. Many superfluid phenomena are 
strongly connected to the phase of the fluid. For instance 
vortices [^j which have been created in a stirring manner 
similar to the rotating bucket experiment |}| , and solitons 
J|,[| , are typical examples of a condensate with a nonuni- 
form phase. Also quantum shock waves causing topolog- 
ical defects in the superfluid, have been created using 
ultra-compressed pulses of slow light M. Here we show 
how the properties of slow light can be used to detect the 
phase of the macroscopic wavefunction of the condensate 
nondestructively, where the slow light receives the phase 
of the condensate in a way similar to phase imprinting 

II- 

The recent experimental demonstrations with slow 
light has revealed a whole new playground with interest- 
ing applications . In this paper we show the advan- 
tages with slow light when the phase of the condensate 
is of interest. Light, and especially the phase contrast 
imaging method, has been successfully used to probe the 
density of the condensate in situ ^ . In order to measure 
the phase, on the other hand, destructive methods have 
been used where the interference pattern has revealed the 
phase of the condensate Jl(J . In this paper we show using 
slow light that the phase of the condensate can be im- 
printed on the light in a nondestructive way, making it 
possible to follow the phase evolution of the condensate 
in situ. 

The envisioned experiment is based on Electromagneti- 
cally Induced Transparency (EIT) |Tl|-[T^] . Suppose that 
the condensate is trapped in an elongated trap and il- 
luminated by a uniform control beam. The beam con- 
trols the group velocity v g of a second beam, the probe 
beam. The probe shall be a pulse with a frequency in 
the laboratory frame that matches exactly the atomic 
transition frequency between the ground state and one 
of the excited states coupled by the control beam. The 
condensate density, which can be monitored in-situ by 
phase-contrast microscopy , is supposed to be station- 



ary, but the phase can in certain cases fluctuate [|l-4| . |l 5[j . 
For low group velocities the probe beam couples strongly 
to the condensate phase gradient which results in a phase 
imprinting on the light. The condensate phase can conse- 
quently be read out by analyzing the interference pattern 
between the initial and final probe light. 

Let us start by deriving the equations of motion for 
the slow light and the Bose-Einstein condensate. The 
probe light is described by a real scalar field tp where we 
ignore the polarization, and ip represents the electric field 
in units of the vacuum noise, E — (h/eo) 1 ^ 2 ^^- Here uiq 
denotes the resonance frequency of EIT. In a medium at 
rest, slow light is subject to the principle of least action 
with the Lagrangian density Q 

Cl = \ ((1 + a)(d t p) 2 ~ c 2 {V V ? - a Wo V) (1) 

where the group index a corresponds to a group velocity 
of 
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and a is proportional to the density of the condensate, p, 
and inver sely proportional to the intensity of the control 
beam Mil, 
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Here the intensities of the probe and control fields are 
calibrated in terms of the Rabi frequencies fl p and f2 c , re- 
spectively. The control beam will dominate and in prac- 
tice \n p \ 2 /\fl c \ 2 does not exceed 1CT 1 0. This means 
the less intense the control beam is the slower the light 
is. The corresponding momentum density of the light 
which is needed when describing the coupling between 
the light and the condensate, is obtained from the sym- 
metric energy-momentum tensor p7| , ^2| 
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Consider for the time being the general situation with 
a moving condensate. The condensate with the flow u 
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and the probe light must now be taken as a combined 
dynamical system with the total Lagrangian density 
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where the condensate is described by the Gross-Pitaevskii 
Lagrangian density 

£M--p(^+yu 2 + £(VVp) 2 + |p + y). (6) 

Here m is the atomic mass, g characterizes the atom- 
atom collisions, V denotes the external potential and p 
the condensate density. The coupling between light and 
matter is given by the relation Ea] 



mu = h\7S + a Q P, 
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where ao is a constant. This coupling gives the correct 
equation of motion 



d tP + V(pu)=Q 
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from the Euler-Lagrange equation. 

Using the Lagrangian in Eq. (^) we derive the equation 
of motion for the light field. The resulting wave equation 
is 

((1 + a)dt - c 2 V 2 + coal + ad t u • V 

+V • aud t )(f{r) = 0. (9) 

In order to simplify things we assume the light field can 
be expressed as a propagating pulse with a single fre- 
quency ujq and a slowly varying amplitude. The resulting 
equation of motion is then of the Schrodinger type 
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where the coupling between light and matter is primarily 
described by the phase gradient of the condensate. 

Let us illustrate the consequences of the coupling be- 
tween the light and the phase of the condensate by using 
an elongated quasicondensate where the phase can fluctu- 
ate JlJJl^l • Fluctuations of the density and the phase of 
a condensate are related to the elementary excitations. 
The density fluctuations are dominated by the excita- 
tions of the order of the chemical potential p. If the con- 
densate is very elongated but still of a 3D character, the 
wavelength of the density excitations are much smaller 
than the radial size of the condensate. The fluctuations 
are therefore of an ordinary 3D form and are small. Con- 
sequently the total field operator can be written as 



*(r) = VPo(r)e 
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where po( r ) is the stationary density and 0(r) is the op- 
erator describing the phase of the condensate. This op- 
erator is given by Q 



0(r) = 
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where a„ is the quasi particle annihilation operator with 
quantum number v and energy e„. The mode functions 
/+ = u v + v u are the sum of the two functions u v and 
v v which are the solutions to the Bogoliubov-deGennes 
equations. 

In an elongated condensate the excitations are mainly 
of two kinds: axial excitations with e v < huj p and ra- 
dial excitations e u > Uu p . The latter has a 3D character 
since the wavelength is typically less than the radial size 
of the cloud and consequently the fluctuations are small. 
The axial excitations on the other hand have wavelengths 
larger than the radial size of the cloud and have a ID 
behavior. Therefore these excitations will be most im- 
portant for the axial fluctuations of the phase. In order 
to actually calculate the phase we note that for a har- 
monic external potential the low energy axial modes |2C[ | 
are described by the energy e v = ^hu z y/ v(y + 3) and 
the functions 



,+ = / Q + 2)(2z/ + 3)gpo(r) (1<1)( z 



where Pu'^ are the Jacobi polynomials, R the radial 
size and L the axial length of the cloud. The quasi- 
particle annihilation operators are now replaced by com- 
plex amplitudes 7 and 7*. To reproduce the quantum 
statistical properties of the phase the amplitudes are 
sampled as random variables with a zero mean value, 
(71/) = (7^) = 0, and bosonic number density 
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where j3 is the inverse temperature fl4[| . 

If the coupling is dominated by the phase gradient, in 
other words we neglect the nonlinear term for the light 
in Eq. (|ltj|), it is clear that we can use the light as a weak 
probe which is affected by the condensate phase gradient. 
The condensate will necessarily be very elongated. We 
can therefore solve the dynamics in ID since the low en- 
ergy axial excitations will acquire a ID character. From 
Eq. ([l0]) it is immediately clear that the light will pick 
up a phase 



ip(r, z) — ip(r, z)e l ^ z ^ 
where the the phase is given by the integral 
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with ko = — and S(z) the phase of the condensate. This 
scenario requires that the interaction between the probe 
light and the condensate is turned on suddenly. This can 
indeed be achieved by tuning the control beam such that 
the probe beam is rapidly slowed down and allowed to 
propagate in the condensate. This results in a pulse delay 
typically of the order of a few micro seconds (t]] . 
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FIG. 1. a) The experimental setup with the control beam 
and the probe beam. The probe beam reveals the phase prop- 
erties of the condensate, b) The probe beam is split where one 
part is let to interact with the condensate. Detection of the 
condensate phase is made with an interference experiment. 
Note the transversal interference pattern. 

To illustrate the phase measuring process we solve Eq. 
( |l0| ) for a light pulse moving with the group velocity 
1mm/ s during 1.5ms in a condensate with N = 3.3 x 10 4 
Na atoms, uj z = 2tt x 10Hz and uj r = 2n x 400-ffz. We 
choose a low group velocity and a relatively long prop- 
agation time in order to show any effects on the light 
intensity from the condensate phase. In Fig.(|^) the con- 
densate together with one possible phase representation 
is shown at a temperature well below the condensate crit- 
ical temperature, T = 0.75Tbec- From Fig. (2a) it is 
clear that only the low excitations are significant, as ex- 
pected. The phase fluctuations are clearly present and 
are typically of the order of 7r/3 in the center of the con- 



densate. The shape of the light pulse, in Fig. (2b), is only 
slightly altered due to the acquired condensate phase. 
The pulse shape changes are typically small compared to 
the overall pulse width and are therefore not necessarily 
very well suited for detecting any acquired phase from 
the condensate. In Fig. (2c) we show the resulting inten- 
sity when the light which has acquired the phase of the 
condensate, and the initial probe beam interferes. Note 
here that the interference pattern will be in the trans- 
verse direction compared to the pulse propagation. For 
sufficiently low group velocities the acquired condensate 
phase is amplified by the factor hko/mv g , which makes 
it possible to study also small phase fluctuations not de- 
tectable by a time of flight measurement H|. As in any 
measurement process, there is inevitably going to be a 
back-action on the condensate from the light. This back- 
action can, however, be made arbitrarily small by choos- 
ing an appropriate ratio |Q p |/|fi c |. The group velocity, 
which is a function of the control beam, is not affected 
by this ratio. 

The ID character of the condensate-light phase im- 
printing has here been used for simplicity. Similar phase 
fluctuations in 2D condensates can also exist |21| and 



consequently the same situation emerges concerning the 
properties of the light. One important factor, however, 
is the inherent time dependence of the phase fluctuations 
in the condensate. As was already seen in Fig. (|^) the 
phase fluctuations are dominated by the lowest modes. 
This means also that the dynamics is governed by the 
corresponding frequencies uj v = e v /%. For a typical ex- 
ternal harmonic potential with lu z = 2ir x 10Hz we ob- 
tain a timescale for the phase fluctuations of the order of 
20ms, corresponding to the ten lowest modes, which is 
long enough to treat the phase as stationary in the im- 
printing process. The typical timescale to apply and slow 
down the light pulse is of the order of micro seconds (tJ 
which allows us to assume an instant phase imprinting. 

In summary we have shown that light with extremely 
low group velocity can be used to probe the phase of a 
condensate. The method was illustrated by imprinting 
the fluctuating phase of a quasicondensate onto the slow 
light propagating through the condensate. The phase im- 
printing technique also allows for studying other forms 
of phase gradients such as vortices [E2j and solitons [E| . 
Especially solitons in two-component condensates with 
repulsive interaction could be detected without opening 
the trap Two-component solitons are significantly 
more difficult to observe because the total density will 
always stay constant, compared to the single component 
condensate soliton which shows a density notch at the po- 
sition of the soliton. The probe light in this case would 
simply pick up the phase of the soliton solution. If the 
light propagates during a long time in a medium with 
a nonuniform phase it is also possible to study scatter- 
ing of the light from the condensate phase [p| which in 
principle can also be used to measure the phase gradi- 
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ent, although this method will not necessarily reveal the 
condensate phase as clearly as a direct interference ex- 
periment. 
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FIG. 2. (a) One possible phase at the temperature 
T = 0.75Tbec ■ Only the lowest modes are important in 
the phase fluctuations, (b) The shape of the light pulse (solid 
black line) is only slightly altered due to the acquired phase. 
The grey curve shows the initial pulse shape (arb. units). The 
dashed line shows the density of the condensate (arb. units) 
(c) The resulting intensity of the light when interfering with 
the initial probe light (arb. units). 
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